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Abstract. The evolution of the line width - luminosity relation for spiral galaxies, the Tully-Fisher relation, strongly constrains 
galaxy formation and evolution models. At this moment, the kinematics of z > 1 spiral galaxies can only be measured using 
rest frame optical emission lines associated with star formation, such as Ho- and [O in] 5007/4959 and [On]3727. This method 
has intrinsic difficulties and uncertainties. Moreover, observations of these lines are challenging for present day telescopes and 
techniques. Here, we present an overview of the intrinsic and observational challenges and some ways way to circumvent them. 
We illustrate our results with the HST/NICMOS grism sample data of z ~ 1.5 starburst galaxies. The number of galaxies we 
can use in the final Tully-Fisher analysis is only three. We find a ~ 2 mag offset from the local rest frame B and R band Tully- 
Fisher relation for this sample. This offset is partially explained by sample selection effects and sample specifics. Uncertainties 
in inclination and extinction and the effects of star formation on the luminosity can be accounted for. The largest remaining 
uncertainty is the line width / rotation curve velocity measurement. We show that high resolution, excellent seeing integral field 
spectroscopy will improve the situation. However, we note that no flat rotation curves have been observed for galaxies with 
z > 1 . This could be due to the described instrumental and observational limitations, but it might also mean that galaxies at z > ! 
have not reached the organised motions of the present day. 
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1. Introduction 

The Tully-Fisher relation (hereafter TFR) is a tight empirical 
relation between the flat rotation curve (RC) velocity and the 
luminosity of spiral galaxies (Tully & Fisher [T9"77l ). The TFR 
has been used as a distance estimator and for measurements of 
the Hubble constant H (e.g. Tully & Pierce I^UUOl . 

In addition to its empirical applications, the TFR is inter- 
esting in itself because it defines a tight relation between the 
total (dark matter dominated) mass of spiral galaxies and their 
luminosity produced by baryons. Furthermore, assuming a stel- 
lar mass-to-light ratio M*/L, the stellar masses of galaxies can 
be calculated from their luminosities and a stellar mass TFR 
can be derived (Bell & De Jong 2001). After addition of the 
gas mass, one obtains the baryonic mass TFR (Verheijen 2001 
Bell & De Jong l^Olfll McGaugh et al. 2000). The variations in 
the dark-to-baryonic mass ratio of galaxies are small and devia- 
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tions from the baryonic TFR are absent down to very low mass 
galaxies (McGaugh et al. 2000 1 although others claim a slight 
deviation for dwarf spirals (Stil & Israel 1998 1. Semi-analytical 
models of galaxy formation struggle to explain simultaneously 
the slope, zero point and tightness of the TFR in all optical and 
near infrared bands (Van den Bosch|3UU3j). 

The tight fundamental relation between mass and luminos- 
ity is interesting to study in the context of galaxy evolution. 
The study of the evolving TFR with redshift can provide valu- 
able information on the luminosity evolution of galaxies and 
the buildup of stellar mass as a function of galaxy mass. As 
we will show in this paper, the analysis of high redshift TFRs 
needs careful treatment of observational limits, selection ef- 
fects, sample definitions and starburst influences, and high res- 
olution high quality spectra. 

In the local universe, HI is used to measure the velocity pro- 
files of spiral galaxies. The gas disk in spiral galaxies extends 
2-3 times further out than the stellar disk. HI measurements 
are currently limited to low redshift, beyond redshift ~ 0.2 
HI emission has not been observed (Zwaan et al. 2001) and 



2 



L. van Starkenburg et al.: On measuring the Tully-Fisher relation at z > 1 



one has to rely on other kinematic tracers. Another gas tracer 
could be CO, but CO detections in the high redshift universe 
are currently limited to very CO bright galaxies, like submil- 
limeter galaxies, quasi stellar objects (QSO) and radio galaxies 
(Hainline et al. 2004 and references therein). One gravitation- 
ally lensed Lyman Break Galaxy (LBG) has been detected in 
CO (Baker et al. 20043. A second attempt to detect a LBG 
in CO failed, although the attempt was on the dustiest LBG 
known (Baker et al. I2004b> . 

Bright optical narrow emission lines like Ha, H/3, 
[Om]5007/4959 and [On]3727 can also be used to trace the 
rotation curve of galaxies. Their presence is limited to the stel- 
lar disk (or more precisely: the star forming disk). In the local 
universe, the agreement between HI and HII measurements of 
RCs is excellent (Courteau 1997). However, these lines shift 
out of the optical regime at redshifts 0.4 - 1 .4. In recent years, 
high resolution near infrared spectrographs like the Infrared 
Spectrometer And Array Camera (ISAAC) at the Very Large 
Telescope (VLT) of the European Southern Observatory (ESO) 
have become available, opening the window out to redshift 2.4 
(for Ha) TFR studies. Examples are Rigopoulou et al. (|2002i, 
who studied massive z ~ 0.6 galaxies, Barden et al. (2003 1 
who found an offset from the local B band TFR of ~ lmag at 
Z ~ 0.9, Lemoine-Busserolle et al. (|2003 1 who used gravita- 
tional lenses to study two galaxies at z ~ 1 .9 and Pettini et al. 
(2001 1, who studied Lyman Break Galaxies at z ~ 3. We used 
ISAAC to study the kinematics of a sample of z ~ 1.5 Ha emit- 
ting galaxies and we present the results in this paper as a case 
study for z > 1 TFR studies. 

Our sample is a subsample of the McCarthy et al. ([1999 > 
HST/NICMOS grism survey sample. McCarthy et al. (JT333} 
surveyed 64 square arc minutes with the slitless NICMOS 
G141 grism and detected 33 emission line objects with vary- 
ing 3cr limiting line fluxes down to 1 x 10~ I7 erg s^'cirT 2 . They 
argue that the detected emission lines are Ha between redshift 
0.75 and 1 .9. The Ha+[N n]6548/6584 complex is not resolved 
due to the low spectral resolution (R ~ 150) of the grism and 
therefore contamination by other emission lines (particularly 
[Om]5007/4959) cannot be excluded and no kinematic infor- 
mation is obtained. This sample is biased to galaxies with large 
Ha equivalent width, EW(Ha), and Ha flux, F(Ha), due to the 
low spectral resolution of the grism. We chose this sample be- 
cause it has clear selection criteria and all sources have known 
Ha fluxes. 

Here, we present observations of 9 objects from the 
McCarthy et al. sample with the ISAAC at the VLT in medium 
resolution mode (R ~ 3000-5000). Our aim was to resolve the 
Ha+[Nn] complex (or the [Ora]5007/4959 doublet) to con- 
firm redshifts, measure accurate line fluxes and linewidths and 
if possible also rotation curves. We use this data to present our 
case study for z > 1 TFRs. 

Hicks et al. (20021 also performed a follow-up of the 
HST/NICMOS grism sample. They observed 14 objects aim- 
ing to detect emission lines, particularly [On] 3727 in the opti- 
cal (R/I band) using LRIS at the Keck telescope. They observed 
in low resolution mode (R ~ 350 - 700 depending on the grat- 
ing used) and therefore did not obtain any kinematic informa- 
tion. They confirmed the redshift from McCarthy et al. (1999) 



for 9 out of 14 objects. They explained the non-confirmations 
by twilight observations or the presence of a nearby bright star 
(emission lines may very well be not bright enough to detect in 
these two cases). In two cases, the [On]3727 line was outside 
the observed wavelength range and other emission lines like for 
example C n]2326, Cm] 1909 and Mgn2800 might have been 
too faint to detect. The fifth non-detection was explained by 
reddening or the emission line detected by McCarthy et al. was 
not Ha but HyS or [O m]5007/4959. In the latter case, no bright 
emission lines are expected in the wavelength range observed. 

Our follow-up is complementary in two ways: we observe 
the objects accessible from the southern hemisphere whereas 
Hicks et al. observed from the northern hemisphere. Only one 
object (J093 1-0449) is in both samples. Second, we observe at 
higher resolution, resolving the emission lines. 

This paper is organised as follows. The first part of the pa- 
per describes the case study data set: the observations of the 
NICMOS grism sample (Section 2), data reduction and analy- 
sis (Section 3), the sample properties (Section 4) and notes on 
individual galaxies (Section 5). The second part of the paper 
discusses high redshift TFRs using the earlier discussed dataset 
as an illustration with a strong focus on pitfalls (Section 6). We 
conclude with a summary and conclusion in the final section 
(Section 7). 

Throughout this paper, we assume Q.m = 0.3, Qa = 0.7 and 
Hq = 70 km s _1 Mpc _1 . All magnitudes in this paper are Vega 
magnitudes. 

2. Sample selection and observations 

The McCarthy et al. (Q599I catalog contains 33 galaxies with 
redshifts between 0.75 and 1.9. We selected all targets from the 
McCarthy et al. (1999 1 sample accessible with the VLT and 
with the line falling in the J or H atmospheric window. We did 
not select on morphology or emission line flux. 

The observations were done in two runs. In the nights of 
February 23, 24 & 25 2002 (ESO program ID 68.A-0243(A)) 
we observed 7 targets in visitor mode using the VLT ISAAC 
long slit spectrograph in medium resolution (MR) mode with 
1" slit under varying atmospheric conditions. In winter 2003 
(ESO program ID 70. A-0304(A)) two targets were observed in 
service mode under excellent seeing and sky conditions (seeing 
< 0/6, clear/photometric) with the 076 slit. Integration times 
varied between 120 to 250 minutes, depending on the atmo- 
spheric conditions and the emission line flux. One target could 
not be acquired, although it was attempted several times. An 
overview of all observations is given in Tabled 

The observational set-up was as follows. The slit was 
aligned along the major axis of the galaxy as determined from 
NICMOS H band images. Where possible without deviating 
more than 10° from the major axis of the galaxy, a bright ref- 
erence star was also included in the slit to make sure the slit 
was on target. To facilitate sky subtraction, total integration 
times were dived in 12 or 15 minutes exposure times, nod- 
ding in ABBA cycles along the slit. After observation of each 
object, a bright nearby standard B star was observed with the 
same instrument setup and the same air mass to allow accurate 
flux calibration. For the 076 slit observations, the B star was 
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Fig. 1. Two dimensional spectra, from the upper left and then clockwise: J0627-6512, J0738+0507a, J1143-8036a and 
J0738+0507b. Note the extended emission in J0627-6512 and J0738+0507b and the positions of the OH sky lines in all spectra. 



Table 1. Observations 



Source ID 


run ID 


slit width 


7hit(s) 


seeing h 


J0627-6512 


68.A-0243(A) 


1" 


7200 


a/ 69 


J0738+0507a 


68.A-0243(A) 


1" 


7200 


a' 83 


J0738+0507b 


70.A-0304(A) 


a/6 


15000 


0760 


J093 1-0449 


68.A-0243(A) 


i" 


7200 


a' 80 


J1056-0337 


68.A-0243(A) 


i" 


14400 


0767 


J1143-8036a a 


68.A-0243(A) 


i" 


7200 


07 78 


J1143-8036b a 


68.A-0243(A) 


i" 


7200 


a/78 


J1143-8036c 


68.A-0243(A) 


i" 


7200 


a/ 91 


J1143-8036d 


70.A-0304(A) 


a' 6 


12000 


a' 60 



a Jl 143-8036a and Jl 143-8036b were observed in the same slit. 

b The seeing was measured on the brightest object in the slit in the 
reduced image. In one case, J1056-0337, there was no bright ob- 
ject in the slit and the seeing was measured on the standard star 
for flux calibration. 



also observed with the 2"slit to calculate the (wavelength de- 
pendent) slit loss correction. Depending on wavelength and slit 
width, the sampling was 0.57 to 0.81 A pix~'. The full width 
half maxima (FWHM) of the sky lines varied between 3.3 and 
4.7 A. 



3. Data reduction and analysis 

We used standard eclipse (Devillard 1997) and IRAF proce- 
dures for data reduction. The available twilight flats were used 
to create a bad pixel map. We removed bad pixels, ghosts and 
cosmic rays in all frames before combining them. Except for 
one object (Jl 143-8036c) we used dome flats for flat fielding. 
An illumination correction to these flat fields was done to re- 
move a small residual gradient in the sky. Residual bias subtrac- 
tion was only necessary for the objects with the highest quality 
data (J0738+0507b and Jl 143-8036d). The spectral tilt was re- 
moved using star traces. The OH lines were used to correct for 
the curvature of spectral lines. If the detected emission line was 
close to one of the edges of the detector, we recalculated this 
correction optimising it for the area around the emission line to 
minimise OH line residuals where they are most relevant. 

Flux calibration was done using bright B stars, observed di- 
rectly after the object. A (wavelength dependent) slit loss cor- 
rection was applied to the 076 slit spectra. The OH lines were 
used for wavelength calibration, < 0.5 A residuals remained 
after a third order fit. 

One or more emission lines are immediately visible in the 
two dimensional spectra in 5 out of 9 cases. In two cases, we 
also detect continuum emission, in one case, we detect con- 
tinuum emission without an emission line (J1056-0337). One 
detection turns out to be a Seyfert 1 (J093 1-0449). The re- 
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duced two dimensional spectra of the detected emission lines 
are shown in Fig.^( exce pt the Seyfert). 

We extracted one dimensional spectra by cutting out a strip 
from the two dimensional spectrum containing all flux of the 
detected emission line, or, if there was no (clear) detection in 
the 2D spectrum, a strip was cut out at the expected position 
of the emission line (using the known distance between the ob- 
ject and the reference star). Extracting the spectrum by trac- 
ing the spectrum was not an option, because we detect weak 
continuum emission in only three sources. The spectra were 
smoothed with a Gaussian with FWHM approximately equal 
to the FWHM of the OH lines in the raw frames and are shown 
in Fig. |2 In the one dimensional spectra, a second or third 
emission line is immediately evident in two cases. The bright- 
est emission line of every detected object was clearly visible in 
the two dimensional spectrum. 

Line fluxes and widths were measured by fitting a Gaussian 
to the detected emission lines using IRAF's 'splot'. Errors were 
estimated by repeated fitting with different parameters. If there 
was severe OH line contamination, we interpolated over the 
OH line to correct for flux losses. Error bars are naturally larger 
in this case. Linewidths were corrected for instrumental broad- 
ening and converted to W20, the width at 20% of maximum flux. 
Line fluxes and widths are given in Tables[2]and[3]respectively. 
Values corrected for OH line contamination are marked by the 
subscript 'OH'. 

We calculated rest frame B or R magnitudes (depending on 
redshift) from the observed F110W (J hereafter) and F160W 
(H hereafter) magnitudes. An H band magnitude was available 
for all targets, J band for a subset only. When J band photom- 
etry was unavailable, we used the average J-H color (equal to 
the median color) of the entire NICMOS grism sample to esti- 
mate the J band magnitude. We calculated the rest frame mag- 
nitudes by interpolating between the J and H fluxes, depending 
on redshift this gives us a rest frame B or R absolute magnitude. 
Where we could not interpolate to obtain a rest frame B or R 
magnitude, we made a rough estimate of this magnitude by us- 
ing the closest flux point available. The errors in the absolute 
magnitudes were calculated as follows: while interpolating be- 
tween the J and H magnitudes, interpolating meaning that the 
redshifted effective wavelength of the B or R band lies in be- 
tween the effective wavelengths of the J and H band, we set the 
error in the measured magnitudes to 0.1, and in the J magni- 
tudes calculated from the average J-H color to 0.4 (=scatter in 
J-H color). We then interpolated the fluxes and uncertainties to 
get the rest frame magnitudes and errors. When the effective 
wavelength of the redshifted B or R band was not between the 
effective wavelengths of the J and H band but inside the wave- 
length range of the J or H band, we set the error to 0.5 mag. If 
it was outside the wavelength range of the J and H band, we set 
the error to 2.0 mag. These numbers are a bit arbitrary, but are 
intended to reflect the increased uncertainties in the magnitudes 
estimates. The apparent infrared magnitudes and correspond- 
ing absolute rest frame magnitudes are given in Table|4] When 
calculating the offsets from the TFR, we use only those points 
where the redshifted B or R filter at least overlapped with the 
observed J or H band. 



Table 4. Apparent magnitudes and absolute magnitudes of tar- 
gets 



Source ID 




J 




H 




M B 






Mr 




J0627-6512 




22.2 




20.4 




-22.3 ± 


0.1 


(- 


-23.1 ± 


0.5) 


J0738+0507a 


< 


18.8 


> 


17.9 


(- 


-22.5 ± 


2 ) 




-23.7 ± 


0.1 


J0738+0507b 


< 


22.9 


> 


22.0 




-20.7 + 


0.1 


(- 


-21.5 ± 


0.5) 


J0738+0507b([Om]) 


< 


22.9 


> 


22.0 




-21.7 ± 


0.1 


(- 


-22.9 ± 


2 ) 


J093 1-0449 




19.7 




19.0 


(- 


-22.2 ± 


2 ) 




-23.2 ± 


0.1 


J1143-8036a 


< 


22.3 


> 


21.4 


(- 


-20.6 ± 


0.5) 




-21.4 ± 


0.1 



a J and H band magnitudes from McCarthy et al. 1 1999 1. The J mag- 
nitudes calculated from the H band magnitude and average J-H 
color are marked by < >. The extrapolated absolute magnitudes 
and their errors (see text) are in parentheses. 



4. Results 

We detect one or more emission lines in 5 out of 9 spectra. 
One dimensional spectra (of the relevant wavelength ranges of 
the original ISAAC spectra) are shown in Figs.|5^-e. In these 
figures, the expected positions of the [Nn]6584/6548 lines (as- 
suming the brightest emission line detected is Ha) are marked. 
If the observed emission line is not Ha, the next most likely 
candidate is [Ora]4959/5007. Also marked is the expected 
position of [Om]4959 (assuming the bright line detected is 
[Om]5007) if the detection of the [Nn] doublet is uncertain. 
To avoid confusion between emission lines and OH line resid- 
uals, the sky spectra are also shown. Other possible identifica- 
tions of the emission lines can be ruled out or are far less likely: 
[On]3727 would put the sources at redshifts larger than 3 (H 
band detection) and would be resolved in a doublet which is 
not observed. The equivalent width of H/3 is in general too low 
to be detected in the McCarthy et al. (1999) survey. We find 4 
Ha emitting galaxies and 1 (likely) [Om]5007/4959 emitting 
galaxy. In Sect. 5, we discuss all galaxies individually. 

In Table |U we list the wavelengths and fluxes from 
McCarthy et al. ( 1999 1. We note that there is a systematic offset 
between the wavelength as found by McCarty et al. ( 1999 1 and 
ours, although all our wavelengths lie within 3cr error bars of 
the NICMOS wavelengths. We checked some of the OH lines 
in the ISAAC spectra and they were correct within a few A. 
We also note that the emission line fluxes are not always in 
agreement. This is probably due to a combination of the low 
resolution of the NICMOS grism and slit losses with ISAAC. 

In Table|3]we list starformation rates (SFRs) and dynamical 
masses for all detected objects. SFRs were calculated using 

SFR^yr- 1 ) = L„ ff (erg s - l )/l.26 x 10 41 (1) 

(Kennicutt et al. (1998) for a Salpeter Initial Mass Function 
(IMF)). 

Dynamical masses were calculated using 
RV 2 

M(R) = — (2) 

where the velocity V = W20/2 and the diameter R = D/2. D is 
the diameter of the galaxy measured as the total extent in the 



L. van Starkenburg et al.: On measuring the Tully-Fisher relation at z > 1 
Table 2. Emission lines detected 



5 



source ID 




line ID 


z 


F a 


pub 
1 OH 




J0627-6512 


1.72616 


Ha 


1.630 


4.0 ± 0.2 






J0627-6512 


1.73167 


[Nn]6584 


1.630 


0.59 ± 0.05 






J0738+0507a 


1.1966 


Ha 


0.824 


2.9 ±0.3 


~ 4.1 


-4.7 


J0738+0507a 


1.2005 


[Nn]6584 


0.824 


2.0 + 0.2 


-2.3 


-3.7 


J0738+0507a 


1.19411 


[Nn]6548 


0.824 


1.6 + 0.2 






J0738+0507b 


1.7354 


Hff or[Om]5007 


1.644 or 2.466 


0.74 ± 0.02 


-0.9 


- 1.0 


J093 1-0449 


1.2973 


Ha + [Nll]6548/6584 


0.977 


21 ±5 






J1143-8036a 


1.53137 


Ha 


1.333 


0.73 ± 0.06 






J1143-8036a 


1.53623 


[Nn]6584 


1.333 


0.16 ±0.03 







a Fluxes in units of 10 16 erg s 'cm 2 . 

b Fqh is an estimate of the emission line flux F had it not been contaminated by one (or more) OH sky lines. 



Table 3. SFRs and masses. 



Source ID 


SFR 

M yr' 


SFRoh 
M©yr 


i 


W 20 
kms~' 




W20 OH 

kms~' 




M R a 
10 10 M o kpc 


J0627-6512 


57 ±3 






344 ± 


11 






4.8 7 


J0738+0507a 


20 ±2 


-28- 


32 


398 ± 


38 


- 479- 


677 


5.1 6 


J0738+0507b b 


10.9 ±0.3 


- 13- 


15 


166 ± 


3 


- 216- 


235 


1.3 8 


J0931-0449 c 


167 ± 40 






5300: 


1 1800 








J1143-8036a 


8.2 ± 0.7 






274 ± 


18 






1.7 4 



a R is half the diameter, measured as the total extent along the slit in the spectrum. 
b SFR and mass were calculated assuming the emission line observed is Ha. 
c J093 1-0449 is a Seyfert 1 galaxy, the SFR is not meaningful. 



Table 5. Comparison with the line fluxes and wavelengths of McCarthy et al. 



Source ID 


AmccQm 


:) /i(jura) 


WA) F* McC 


p* 


OH 


J0627-6512 


1.742 


1.72616 


158 


1.8 ±0.5 


4.0 ± 0.2 




J0738+0507a 


1.210 


1.1966 


134 


16 ± 1.5 


2.9 ±0.3 


- 4.1-4.7 


J0738+0507b 


1.77 


1.7354 


346 


0.9 ±0.3 


0.74 ± 0.02 


-0.9-1.0 


J093 1-0449 


1.299 


1.2973 


17 


24 ± 1.7 


21 ±5 




J1143-8036a 


1.538 


1.53137 


66 


1.2 ±0.4 


0.73 ± 0.06 





a All fluxes are in units 10 16 erg s 'cm 2 . 



spectrum. We also measured the diameters in the images with 
gave consistent results. The masses and radii are also listed in 
Tabled Note that these masses are lower limits as no correction 
for inclination or OH lines has been applied. 

The [Nul/Ha ratio can be used to get an estimate of the 
metallicity of galaxies. We used the calibration of Denicolo et 
al. (|5002 1 and the results are reported in Table|6] 

5. Notes on individual objects 

We will now discuss all galaxies individually, paying attention 
to the identification of the emission line(s), Ha/[Nn] ratios, 



linewidths and kinematics. Where we do not detect any emis- 
sion line, we will attempt to give an explanation. 



J0627-6512 A single bright emission line is visible between 
two bright OH lines (see Figs.[T]i andEJi). Although [N n]6584 
emission is not visible by eye in the two dimensional spectrum, 
it is quite obvious in the one dimensional spectrum. Hence, 
we confirm the redshift to be 1.630. The [Nii]6584/Ha' ratio 
is about 0.13, confirming that we are looking at a star forming 
galaxy (Brinchmann et al. 2003, Gallego et. 1997). We do not 
detect continuum emission in the spectrum. 
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Fig. 2. ID spectra for all sources. Assuming the brightest emission line detected is Ha, the expected position of the 
[N n] 65 84/6548 lines are marked. If [Nn]6584/6548 is not or marginally detected, we also marked the expected position of 
[Om]4959 assuming the brightest emission line is [Oiii]5007. Note that [Ora]4959 and [Nn]6548 fall outside the wavelength 
range observed for object J0738+0507b. 



Table 6. Metallicity 



Source ID 



[Nil] /Ha 



log([Nll]/Ha) 12 + log(0/H) 



10627-6512 0.125 ±0.020 -0.90 ±0.16 

I0738+0507a a 0.69 ±0.14 -0.16 ±0.20 

Sl±0.1 <0±0.1 

I0738+0507b < 0.005 b ± 0.005 S -2.3 ± 1 



8.46 ±0.16 
9.00 ±0.16 
5= 9.12 ±0.09 
5= 7.4 ±0.8 



I1143-8036a 0.22 ± 0.06 



-0.66 ± 0.27 8.64 ± 0.22 



a The OH line corrected values are on the second row. Note this 
galaxy is probably a narrow line AGN and the line ratio cannot be 
interpreted as a metallicity effect. 

b This upper limit is based on the bright part of the emission line, 
measured in 2D image (to get best constraint) 



In the NICMOS image, this object looks like an asymmet- 
ric edge-on galaxy. Indeed, we do detect some extended emis- 
sion in the spectrum on the same side of the galaxy as in the 
image. Remarkably, the emission line is not tilted, and there is 
no sign of ordered rotation. The elongated appearance of the 
galaxy could be intrinsic, and not due to an edge-on orienta- 
tion. We might also miss a tilt in the emission line due to the 
observation conditions, see Sect. 16.31 for a discussion of this 
possibility. The optical seeing during these observation varied 
between 0765 and 1724 (median seeing 0775). With current 



data, we cannot distinguish between the two possibilities. We 
would need excellent seeing, better S /N (the S /N of the ex- 
tended emission is poor) data with a smaller slit to determine 
the nature of J0627-6512. 



J0738+0507a This is a very bright source, both in emission 
line and in continuum flux. We detect Ha, [Nn]6548 and 
[Nn]6584 and continuum emission at redshift 0.823. Because 
of its brightness and its compact morphology, it has been sug- 
gested that J0738+0507a hosts an Active Galactic Nucleus 
(AGN) (McCarthy et al. [T999I Hicks et al. 12002b . However, 
our detection of narrow emission lines rules out the possibility 
of a Seyfert 1. Star forming galaxies and AGNs can be dis- 
tinguished from their emission line ratios due to their different 
excitation properties. The most suited line ratio diagram to sep- 
arate star forming galaxies from AGNs is the line ratio diagram 
with log([N n6584]/Ha) on one axis and log([0 iii]/Hj6) on the 
other axis (e.g. Brinchmann et al. 2004). We do not have mea- 
surements of all these lines, but the ratio log([N n]/Ha) can 
identify some (but not all) AGNs. According to Brinchmann et 
al. 0004), all galaxies with log([Nn]/Ha) > -0.2 are AGNs. 
The measured log([Nn]/Ha) ratio for J0738+0507a is quite 
uncertain, because both lines are contaminated by OH line 
emission. The measured value (see Table |6j) is -0.16 ± 0.20. 
The true value is most probably larger (less negative), because 
the [Nil] line is more contaminated then the Ha line. We sug- 
gest that J0738+0507a is likely a narrow line AGN. 
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Fig. 3. Rotation curve of J0738+0507b. The shaded areas cor- 
respond to OH lines. 



J0738+0507b This galaxy was observed under excellent at- 
mospheric conditions with the 076 slit and a total integration 
time of just over 4 hours. McCarthy et al. (1999) refer to it as a 
"putative emission line". We detect a beautifully tilted emission 
line, extending ~ 178 (1 "corresponds to ~ 8 kpc at the source 
redshift, see below) and rising continuously without flattening 
off, see Fig. and Fig. [3] The rotation curve velocity 2V is 
at least ~ 110 km s (total visible extent of the emission line 
without correction for inclination or the OH line cut off (see 
below)). 

The wavelengths of the OH sky lines turn out to be very 
unfortunate: one bright OH line cuts off the detection on the 
short wavelength side, another falls on the middle of the emis- 
sion line, making it hard to put strong limits on emission line 
flux, extent and velocity. Strictly speaking we can only measure 
lower limits. 

What is clear however, is that we do not observe a double 
horned profile. What we observe is a bright part and a much 
fainter part on the long wavelength side. It is possible that a 
similar fainter outer part is also present on the short wavelength 
side, but this is impossible to detect due the the presence of 
the bright OH line at that side of the emission line. Comparing 
the flux as a function of position in the spectrum to the flux in 
the image, and assuming that the equivalent width of the emis- 
sion line does not vary with position, we checked that there is 
some continuum flux in the image at the undetectable position 
in the spectrum. We could therefore be looking at a centrally 
star bursting system, with lower levels of star formation in the 
outer parts. 

As can be seen in Fig. |2j;, there is no sign of [Nn]6584 
and there is no bright OH line close to the expected position of 
[N n]6584. We can rule out an Ha/[N n]6584 ratio smaller than 
190 at the 3cr level (for the brightest part of the emission line, 



assuming constant Ha/[Nn]6584 for the whole galaxy). The 
highest Ha/[Nn]6584 ratios observed for local starburst galax- 
ies are ~ 20 for Blue Compact Dwarfs with some outliers with 
~ 100 (Gallego et al.1997, Brinchmann et al. 2003), therefore 
identification of the emission line as Ha seems unlikely. 

Unfortunately, the detected emission line is close to the 
edge of the detector. Both [Nn]6548 (assuming Ha) and 
[Om]4959 (assuming [Om]5007) fall off the detector, pro- 
hibiting confirmation of [Om]5007/4959. Although McCarthy 
et al. (1999) could in principle have resolved the [Om] dou- 
blet (the separation between the lines at this redshift would be 
about twice their resolving limit), the small line flux immedi- 
ately explains why they could not in this case. We conclude 
that [Om]5007 is the most likely candidate, putting the red- 
shift at 2.466 (instead of 1.644) and making this object one of 
the highest redshift objects with ordered rotation. We note that 
high redshift [O m] emitting galaxies have been misidentified 
as Ha before: Moorwood et al. ( 2003 1 find that most of their 
presumed Ha emitters are [Om]5007/4959 emitters. 



J0931 -0449 This galaxy is identified as an AGN: it has a very 
broad emission line (FWHM ~ 2500 km s~'), confirmed to be 
Ha (and [Nn]6548 and [Nn]6584) by Hicks et al. d2002b who 
detected [On] 3727. 



J1 056-0337 The emission line of this galaxy was not de- 
tected, although the emission line flux reported by McCarthy 
et al. d 1999b . 4.6 x 10~ 16 erg s^'cirT 2 , is by far not the faintest 
in the sample and we do detect continuum emission. Possible 
explanations are sky line contamination (parts of the spectra 
are very crowded with OH lines), extended emission, or a spu- 
rious source in the NICMOS sample (J1056-0337 is the lowest 
redshift object from the NICMOS sample, the emission line is 
near the edge of the wavelength range covered by the NICMOS 
grism, the signal-to-noise of the detection is only 3). This ob- 
ject has also been observed as part of the FIRES survey of 
MS 1054-03 (Forster-Schreiber et al. l2005> . and has photomet- 
ric redshift of 0.4 (Franx, private communication) whereas the 
redshift according to McCarthy et al. is 0.72. 



J1 1 43-8036a This source has the lowest S/N detection in our 
sample. We see a tentative detection of [N n]6584, although we 
cannot rule out other possibilities. Continuum emission is not 
detected. 



J1143-8036b J1143-8036a and J1143-8036b were detected 
by McCarthy et al. ( 1999) as a pair with nearly identical red- 
shifts. Jl 143-8036b was the fainter one of the pair, and here we 
barely detect J1143-8036a (the distance to the reference star 
of the detected emission line ruled out the other interpretation: 
non-detection of Jl 143-8036a, detection of Jl 143-8036b). As 
J1143-8036b has lower emission line flux than J1143-8036a, 
and J1143-8036a was barely detected, J1143-8036b is proba- 
bly below the detection limit. 
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J 1 1 43-8036C This galaxy was not detected. The emission line 
may fall below the detection limit or the line may fall on top of 
a bright OH line. 

J1143-8036d This galaxy was not detected in the spectrum 
despite excellent observing conditions and long integration 
times. We cannot explain this. 

J1120+2323a This galaxy was scheduled in service mode 
under excellent conditions, but could not be identified on the 
acquisition image, is spite of repeated attempts. 

We confirm line emission in 5 out of 9 galaxies, including 
one Seyfert 1, one possible Seyfert 2 and three starburst galax- 
ies. In one case, the emission line is most likely [Oin]5007 
instead of Ha. The results in this paper are therefore based on 
a very small number of sources (three), one of them having an 
uncertain redshift. 



6. The Tully-Fisher relation at z > 1 

We will now turn to the discussion of high redshift TFRs, us- 
ing the NICMOS galaxies as a case study. We will first use this 
data set to present the z ~ 1 .5 starburst TFR without any cor- 
rections whatsoever. Then, we discuss extinction and inclina- 
tion corrections for high redshift galaxies. We will then explain 
our choices for the velocity parameter, the luminosity parame- 
ter and the local comparison sample, and how the results would 
change if other choices had been made. Finally selection effects 
in velocity, magnitude, sample specifics and star formation are 
discussed. The whole discussion is strongly focused on the pit- 
falls of high redshift TFR analysis in order to assess what can 
be attributed to an evolving TFR and what to peculiarities of 
this and other samples. 

6.1. Thez ~ 1.5 starburst TFR 

For the case study, we study the rest frame B and R band 
starburst TFR using W20 on the velocity axis. As a local ref- 
erence sample, we take the Verheijen (2001) B and R band 
TFR with again W20 on the velocity axis. This sample has been 
corrected for inclination and extinction (extinction correction 
recipe from Tully et al. 1998 1. The Verheijen sample consists 
of spiral galaxies form the Ursa Major cluster and is among the 
tightest TFRs in the literature. The velocities are HI measure- 
ments. 

In Figs. |4k and b we plot the local B and R band TFR of 
Verheijen (2001 1 with our redshift ~ 1.5 objects, which were 
not corrected for inclination or extinction. All points lie signif- 
icantly above the local TFR. Without any corrections (and ex- 
cluding the possible narrow line AGN J0738+0507a), we esti- 
mate that the z — 1 .5 TFR lies ~ 2.0 magnitudes above the local 
TFR in B and ~ 1.8 magnitudes in R at log(W2o[km s~']) ~ 2.5. 
For comparison, using ISAAC in a similar way, Barden et al. 
( 2003 1 found an offset of around 1 magnitude at z ~ 1 . 



6.2. Effects of extinction and inclination 

Contrary to the local reference sample used (and most local 
TFRs), our data has not been corrected for extinction and incli- 
nation. We did not attempt to correct our data points for extinc- 
tion or inclination because both corrections are very uncertain 
for our data. Not correcting for extinction means that we mea- 
sure a lower limit to the offset from the TFR. 

Correcting velocities (line widths) for inclination would de- 
crease the measured luminosity offset. To bring the sample on 
the local TFR, the required shift in log(W2o) is 0.29 in B and 
0.25 in R. Assuming random inclinations, the average shift in 
log(W2o) would be 0.3. Therefore, the lower limit to the lumi- 
nosity offset is ~ 0. 

However, the effects of extinction and inclination are not in- 
dependent: more inclined galaxies are more heavily extincted. 
Depending on inclination and velocity width (if the extinction 
correction is velocity dependent like the Tully et al. ( 1998i) cor- 
rection), the slope of the shift in position between the location 
of the uncorrected and corrected point in the TF plot, is steeper 
or shallower than the slope of the (local) TFR, and the mea- 
sured luminosity offset increases or decreases respectively. For 
large inclinations (edge-on galaxies), the net effect is a lower 
limit on the luminosity offset. In Fig. [5] we plot the combined 
effect of extinction and inclination correction using the (local) 
correction recipe of Tully et al. (1998 1. Plotted are the local 
TFR from Verheijen (2001 ) and vectors showing where points 
on the TFR would lie if they had not been corrected for ex- 
tinction and inclination for inclinations of 80, 70, 60, 50 and 
40°. The figure shows that for galaxies at the relevant velocity 
range, we measure a lower limit to the luminosity offset if they 
are more inclined than ~ 60°. Increasing the amount of extinc- 
tion (and/or a shallower local TFR, the Verheijen TFR is among 
the steepest known) lowers this turnover inclination. One mag- 
nitude extinction more than locally lowers the turnover incli- 
nation to ~ 40°. For a sample of galaxies with random inclina- 
tions, the net effect on the measured luminosity offset is then 
~ 0. Although we cannot measure the inclinations accurately, 
we can say that J0627-6512, J0738+0507b and J1143-8036a 
are significantly inclined, given their elongated morphology 
even at limited spatial resolution. Therefore, we are confident 
that the measured offsets are lower limits to the luminosity off- 
set, although the precise amount remains uncertain due to the 
small number of sources. 

Of course this result relies on the assumption that local ex- 
tinction corrections apply for the galaxies in our sample, which 
is uncertain. Under this assumption however, the measured off- 
set from the local TFR is robust against the combined effect of 
inclination and extinction. 

Accurate inclinations for high redshift galaxies can be ob- 
tained with high resolution imaging. Without any extinction 
correction, the measured offsets from the local TFR are then 
lower limits by definition. If high resolution imaging is not 
available, the above approach minimises uncertainties. For all 
high redshift TFR samples, one has to assume local extinction 
corrections apply for high redshift galaxies unless an indepen- 
dent estimate of the extinction is available. The Balmer decre- 
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Fig. 4. The local B band (left figure) and R band (right figure) TFR of Verheijen (2001 1 with over plotted the z ~ 1.5 NICMOS 
galaxies. The two solid symbols both correspond to J0738+0507b for different redshift identifications. They are connected by an 
arrow to make this clear. The triangle is J0738+0507b assuming the detection emission line is [Om] instead of Ho- (this point 
is not in the R band TF plot because its rest frame R magnitude is highly uncertain, see Table@J. J0738+0507a is presumably a 
Seyfert 2 galaxy and is marked by a cross. The arrows indicate an estimate of underestimation of the linewidth due to influence 
of nearby OH lines. As discussed in Sect. 16.31 the linewidths are subject to a large number of uncertainties. For example, if a 
correction for beam smearing could be made, all linewidths would shift even further to the right. If the error bars were comparable 
or smaller than the point sizes, they were left out for reasons of clarity. 



ment or SED fitting could constrain the extinction, although the 
latter method suffers from degeneracy between age and dust. 

6.3. Velocity: W 20 

A number of different parameters may be plotted on the ve- 
locity axis of the TF-plot. Examples for rotation curves are the 
maximum rotation curve velocity V max and the velocity of the 
flat outer part of the rotation curve VW When using linewidths 
(double horned profiles) one usually takes the width at some 
percentage (50% or 20%) of the maximum flux or average flux 
of both horns. A number of possible corrections may be applied 
to those linewidths, for example for instrumental broadening 
and mass motions. All TFR velocity parameters correlate well 
with eachother. The tightest TFRs are found using Vfl at , which 
traces best the total mass of the galaxy, on the velocity axis 
(Verheiien EPnTl . 

We use the width at 20% of maximum flux W20 as an ap- 
proximation for twice the flat rotation curve velocity Vfl at . We 
correct the measured linewidth for instrumental broadening. 
We choose W20 as our line width definition because in local 
Tully-Fisher samples, the difference between 2Vfl at and W20 is 
smallest (compared with other line width definitions). 

Although there are some prescriptions to convert the HII 
W20 linewidth to velocity or HI linewidth (see Rix et al. 1 19971 
Pisano et al. 2001 1, we did not apply any of those as they give 
contradicting corrections factors. Pisano et al. (2001 1 show for 
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Fig. 5. The combined effect of extinction and inclination. 
Plotted are the local TFR of Verheijen (2001) with the vec- 
tors showing where a galaxy lying on the TFR would lie if it 
had not been corrected for inclination and extinction (accord- 
ing to Tully et al. [T998l for inclinations 80°, 70°, 60°, 50°and 
40°from bottom to top. 
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a sample of nearby blue compact galaxies that the correction 
factor from W 20 (HII) to W 20 (HI) is large (~ 20% or more) 
for galaxies with small linewidths (W2o(HII) £ 140 km/s) and 
small for galaxies with large linewidts (less then 10% for galax- 
ies with linewidths comparable to ours). On the other hand, 
they also find that for galaxies with equivalent width compa- 
rable to ours, the HII linewith is up to a factor 2 smaller than 
the HI linewidth. Rix et al. (1997s cite different studies inves- 
tigating the relation between optical linewidths and flat circu- 
lar velocity. They conclude that l¥2o(HH) should be corrected 
upward by 14% regardless of lV2o(HII), equivalent to reduc- 
ing the luminosity offset by 0.2 mag in B band for the slope 
of the Verheijen (2001 1 TFR. We conclude that the uncertainty 
in the correction factor is significant although the correction 
in log(W2o) (and hence luminosity offset) is probably small. 
Therefore, we decide not to correct the W20 linewidth (except 
for instrumental broadening). 

A crucial assumption in our analysis is that the veloc- 
ity widths and the rotation curve velocity (in the case of 
J0738+0507b) trace the mass of the galaxy in a similar way 
as in the nearby universe. That is, that linewidth traces the or- 
dered rotation of the outer parts of the galaxy where the rota- 
tion velocity is constant. This assumption is highly non-trivial 
to prove (or disprove). Both observational effects and funda- 
mental questions underly this assumption. 

Inflows, outflows, mergers and other dynamical distur- 
bances may bias the line width low or high. Slit effects may 
prohibit the identification of peculiar kinematics. More uncer- 
tainties arise due to slit spectroscopy: the kinematic and photo- 
metric major axes need not coincide, and the velocity profile is 
smeared out because the slit and the galaxy are of comparable 
size. Poor spatial resolution can bias the RC velocity low up to 
a factor two (Erb et al. 2004 1. Poor spectral resolution causes 
smearing in the spectral direction. Poor signal-to-noise, S IN, 
observations cause a significant uncertainty in the linewidths 
and RC velocity observed (Moorwood et al. 2003 1. The com- 
bined effect on the measured line widths is uncertain and can- 
not be quantified with our current data, leaving significant un- 
certainties in the measured offset from the TFR. 

The distribution of the emission line flux over the galaxy 
also influences the linewidths and RC velocity measurements. 
Barton & Van Zee ( 2003 1 showed that a central (or more gen- 
eral, a local) peak in star formation may bias the linewidth mea- 
sured low by a factor 2. 

A more fundamental question is: do the optical emission 
lines extend out to the flat part of the RC like they do in the local 
universe? Although tilted emission lines have been observed 
out to redshift 3.2 (Moorwood et al. 120031 Pettini et al. l2001i . 
no flat RCs have been observed above redshift ~ 1 (Moorwood 
et al. 120031 Pettini et al. l200Tl Lemoine-Busserolle et al. 120031 
van Dokkum & Stanford l200Tl Erb et al. l2003ll2004l . A large 
sample of z > 1, high S/N, high spectral resolution and ex- 
cellent seeing observations may show these disks do exist at 
these redshifts, but it may also be the case that star forming 
disks have undergone significant evolution since z ~ 1 and do 
not extend to the flat part of the RC. Perez ( 2004 1 recently re- 
ported evidence for stellar disk truncation in the redshift range 
0.6 <z < 1.0. 



Kannappan et al. ( 2005 1 investigated the importance of var- 
ious effects on the offset from the local TFR of two high red- 
shift TFR studies (z ~ 0.34 and ~ 0.52). When all samples 
are converted to the same cosmology and reference sample, 
the largest correction factor is that of rotation curve trunca- 
tion (0.7 1 magnitude using their cosmology and reference sam- 
ple). Their samples used RCs only, the situation for linewidths 
might be different, nonetheless this should be a major warning 
that our velocity measurements cause a large uncertainty in our 
analysis. 

6.4. Rest frame magnitude 

TFRs can be measured in all optical and near infrared bands. 
The scatter in the TFR decreases with longer wavelength be- 
cause of decreased sensitivity to dust and star formation (e.g. 
Verheijen 2001 1. Our choice for rest frame B and R band was 
forced by the available photometry. Spitzer will open the win- 
dow to rest frame K band photometry. 

6.5. Local comparison sample 

One can only measure a luminosity offset with respect to a ref- 
erence sample. The galaxies in the reference sample may be 
and probably will be different from the high redshift galaxies 
in many respects. They may have different SFRs and star for- 
mation histories (SFH), higher ages, different metallicities and 
dust properties, different emission line fluxes and EWs, more or 
less or other kinematic disturbances, lower gas mass fractions, 
and more. All these factors contribute directly or indirectly to 
the position of the galaxies in the TF plot and the luminos- 
ity offset measured. Careful definition of the reference sample 
is crucial for the interpretation for the results. In our case, we 
choose to take the Verheijen ( 2001 1 sample as a reference sam- 
ple. This sample contains only 'very ordinary' cluster spirals. 
Strictly speaking, we have measured the offset of young star 
forming galaxies from the 'most ordinary local spirals' TFR. 
We will show now that our conclusions do not change if we 
choose an other reference sample. 

A local starburst sample might be an obvious choice for a 
reference sample. Several local starburst samples for TFR ap- 
plications are available and they follow the normal local TFR, 
with outliers and increased scatter mainly at the low mass end 
(Mendes de Oliviera et al. 120031 Barton et al. 1200 II Van Driel, 
Van den Broek & Baan IT9951 Brungardt IT9881 Davoust & 
Contini 2004 1. The effects of star formation on the TFR be- 
come apparent at log(W2o) ^ 2.4. When studying high mass 
galaxies (like we do), taking a starburst TFR does not make a 
difference for the measured luminosity offsets. We discuss the 
effects of star formation on the TFR extensively in subsection 
IO 

We might also choose a sample containing all local spi- 
rals, following Kannappan et al. (2005 1 who used a morphol- 
ogy blind sample of bright emission line galaxies brighter than 
Mb = -18 and inclined more than 40 degrees. This kind of 
TFR generally has a shallower slope than other TFRs, because 
of the larger effect (see subsection !6.7> of a starburst on lower 
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mass galaxies. The approach also has it pitfalls, as the high red- 
shift sample does not contain 'all' high redshift galaxies, and 
depending on the high redshift sample, the amount of luminos- 
ity evolution may appear different. Using this TFR to measure 
the offset from the local TFR at log(W) = log(2V) = 2.5 re- 
duces the measured offset by 0.2 mag (Kannappan et al. rela- 
tion) plus 0. 15 mag (to account for the different value used for 
Ho) = 0.35 mag. 

Our sample is selected on EW(Ha) and F(Ha), a local sam- 
ple selected on these properties could also serve a reference 
sample. Unfortunately, there is no such sample available. 

We conclude that the differences between local TFRs for 
different samples of galaxies are important for the low mass end 
of the TFR. The differences for high mass end of the TFR are 
negligible compared to the other uncertainties in the analysis. 

6.6. Selection effects 

Selection effects play an important role in high redshift TFR 
studies. We will now discuss selection in velocity, selection in 
magnitude, and selection in equivalent width. The latter is spe- 
cific for our sample. 

In Fig.|S] we plot once and twice the for formal slit limits 
for ISAAC velocity resolution (1" slit, MR mode R ~ 3000) for 
Ha at redshift 1 .5. These limits vary very little between redshift 
1 and 2.5 for Ha measurements with ISAAC. These limits are 
upper limits to the actual resolution because seeing is the limit- 
ing factor. It is immediately evident why measuring TF slopes 
is so difficult at high z: the range in velocity is small compared 
with local samples, where velocity resolution limits only play 
a role in samples of dwarfs. Therefore, we did not attempt to 
measure a 'relation' but only an offset from the local TFR at 
W20 ~ 320 km s corresponding to log(W2o[km s ]) ~ 2.5. 

It becomes increasingly difficult to measure TFR parame- 
ters for galaxies of a certain luminosity with increasing red- 
shift, pushing one again to larger mass (and hence brighter) 
galaxies. There is no clear absolute magnitude cut in the 
NICMOS grism sample, it is biased against low equivalent 
width emission lines (see below). Moreover, the redshift range 
is quite large and the depth of the observations varies from field 
to field. To give an impression of a typical luminosity limit, the 
magnitude limit for the lowest luminosity object in our sam- 
ple (in the total NICMOS grism sample there are a few slightly 
fainter objects) at redshift 1.5 would be -20.4 in B. At the ve- 
locities observed, we cannot distinguish between a zero point 
offset or increased scatter in the TFR in magnitude limited sam- 
ples. 

The main selection effect in the NICMOS grism sample is 
a bias against emission lines with low equivalent width. A se- 
lection effect on EW cannot be translated to a simple straight 
line in the TF plot, like the selection effect on velocity and 
magnitude. However, the selection effect on EW could be very 
important. Kannappan et al. (2002i showed that in the local 
universe EW(Ha) correlates with offset from the TFR, with 
large EW(Ha) galaxies being up to 2 mag brighter in R than 
expected from the local TFR although with considerable scat- 
ter. Kannappan & Barton ( 2004 1 showed a similar result for the 
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Fig. 6. Selection in velocity for the high redshift TFR. Plotted 
are the local TFR of Verheijen ( 2001 1, the data points presented 
here (using the same symbols as in Fig.|4j and data from var- 
ious other authors whose data were also (partially) obtained 
with ISAAC. Filled symbols are RC velocities, open symbols 
linewidths and point size increases with redshift. The vertical 
lines are once (thick line) and twice (dashed line) the ISAAC 
velocity resolution ( 1 " slit MR mode) for Ha at redshift 1.5. 

B band residuals. However, there are three reasons why these 
results should not be simply copied to our galaxies. 

First, the galaxies in those samples have rest frame 
EW(Ha) £ 60 A, and most have EW~ 20 A. The rest frame 
EW(Ha) of our sample of galaxies ranges from 161 to 310 
A (detections only, and excluding the Seyfert 1) (numbers from 
McCarthy et al. U999> . i.e. an order of magnitude larger. A local 
sample of galaxies with comparable EW(Ha) and L(Ha) line 
strengths is unavailable due to the rarity of these galaxies in the 
nearby universe (James et al. 2004 Gallego et al. !1997l >. Based 
on the EW(Ha) alone, we might expect the NICMOS galaxies 
to be brighter than the not EW selected z ~ 1 .5 TFR although 
the amount of brightening remains uncertain. 

Second, Kannappan et al. (2002) and Kannappan & Barton 
(2004) also showed that B-R color correlates with residuals 
from the B and R band TFR, with the residuals being about 
~0.5 mag at the estimated B-R color for the NICMOS galaxies. 
This amount of brightening is significantly less than expected 
from the EW(Ha), indicating again we are looking at an 'in- 
comparable sample'. 

Third, the systematic offset from the TFR correlating with 
EW(Ha) and B-R color as found Kannappan & Barton (2004) 
mainly comes from galaxies with log(W[km s~']) 55 2.4 (their 
definition of line width W is different from ours, but the dif- 
ferences between different definitions of W are minor). The 
galaxies in Kannappan et al. (2002) cannot easily be traced 
through their diagrams. EW(Ha) may be important for low 
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mass galaxies, but irrelevant for high mass galaxies, and hence 
the NICMOS galaxies. 

We therefore conclude that although the selection of high 
EW galaxies most likely biases our result, the effect cannot be 
quantified by a comparison to local samples. We therefore take 
a different approach in the next section. 

6. 7. The effect of a starburst on luminosity in the TF 
plot 

The study of the high redshift TFR is limited to galaxies with 
a sufficient SFR to detect and resolve emission lines like Ha, 
[Om]5007/4959 and [On]3727. These galaxies may not be 
representative (from a TFR point of view) for galaxies with 
smaller SFR. Star formation may affect the position of a galaxy 
in the TF plot in both luminosity and velocity. The effects on 
velocity were discussed in Sect. 16.31 A starburst increases the 
luminosity of a galaxy, and one might worry that the use of only 
starburst galaxies will mimic luminosity evolution in the TFR. 
We will now show that the increase in log(L) (the parameter in 
the TF plot) is negligible for massive galaxies assuming they 
have a significant older population, which - as we will show 
below - might not be the case for some of the NICMOS galax- 
ies. 

We already discussed a hint that star formation could be 
important in Sect. 16.61 namely the results of Kannappan et al. 
( 20021 120051 and Kannappan & Barton (2004i. These papers 
showed that offsets from the U, B and R band TFR corre- 
late with the star formation indicators B-R color and global 
EW(Ho'), bluer color and larger EW (hence more actively star 
forming) galaxies being overluminous for their RC velocity. 
We also noted that different parameters gave conflicting results 
for our sample, and outliers mainly occur at the low mass end 
of the TFR. On the other hand, there was also a hint that star 
formation may not be important for TFR studies. As we saw 
in Sect. 16.51 local starburst TFR follow the local normal spiral 
TFR, with some outliers at the low mass end. We will now try 
to explain these observations using a simple model. 

We plot in Fig. the effect on luminosity of local galax- 
ies for different starbursts (different SFRs, different duration of 
the starburst) calculated using models from Bruzual & Chariot 
(2003 (BQ) (Padova 1994 models, Cabrier IMF, solar metal- 
licity). We added only the luminosity of the starburst, and ne- 
glected the effect of the additional mass on the RC velocity. As 
can be seen in this figure, the luminosity of low mass galaxies 
is dominated by even modest amounts of additional star forma- 
tion. For galaxies with log(W2o) ~ 2.4 only strong starbursts 
have a significant effect on log(L). This is the W20 value where 
the BC model (see Fig.|7J for a short (10 7 yr) 10 M Q yr~' star- 
burst lies 1 mag above the (local) TFR. This figure immediately 
explains why outliers are often low mass galaxies and why the 
zero point of the TFR at the high mass end is so stable. 

High redshift galaxies have a younger and less massive 
population compared with local galaxies. The net effect on lu- 
minosity using BC models with an exponentially decreasing 
SFR is an increase in Lb for high redshift galaxies. Additional 
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Fig. 7. The luminosity of starbursts of different SFRs (indicated 
in the plot) and different durations (10 7 , 10 8 and 10 9 yr for the 
dotted, dashed and long dashed lines respectively) added to the 
local Verheijen (2001 1 TFR (thick line). The thin line indicates 
one mag brightening. 

starformation will therefore have a smaller effect on the lumi- 
nosities of high redshift galaxies compared to local galaxies. 

We conclude that provided the older stellar population is 
sufficiently luminous and the measured velocity represents the 
mass of the galaxy, a starburst will not change the position of 
a galaxy in the TF plot and starburst galaxies are (from a TFR 
point of view) representative for the total galaxy population. 

We will now attempt to check if the NICMOS galaxies ful- 
fill these criteria and are therefore TFR-representative for the 
entire z ~ 1.5 spiral galaxy population. The linewidths in- 
dicate that the NICMOS galaxies are sufficiently massive for 
TFR studies: their linewidths (without inclination correction) 
are log(W) t 2.4 . The non-detection of rotation curves can 
be explained by poor seeing. Improving the velocity measure- 
ments remains the biggest challenge for all high redshift Tully- 
Fisher studies. 

For the case of the NICMOS galaxies, the luminosity of 
the older population is hard to constrain. The SFRs indicate we 
cannot assume the starburst luminosity to be negligible. The 
large EWs indicate the importance of a young stellar popula- 
tion. We used the starburst99 models (Leitherer et al. 119991 to 
estimate the age of the galaxies. Assuming constant star for- 
mation and a Salpeter IMF with mass range 1 - 100MO, the 
EW indicate that the age of the galaxies is ^ 1 Gyr. This results 
strongly depends on the initial mass function and star formation 
history. A younger starburst could also produce the large EW, 
while contributing less to the total luminosity. Without addi- 
tional data, we cannot calculate a reliable estimate of the star- 
burst luminosity, and we cannot answer the question whether 
this sample is TFR-representative for the entire z ~ 1 .5 spiral 
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galaxy population. Note that this conclusion is specific for this 
sample and not for other samples at similar redshift. The EWs 
of this sample are extreme, and we have no photometric data 
that could provide constraints on the SFH. 

7. Summary and conclusion 

We studied the challenges in measuring z > 1 TFR using emis- 
sion line galaxies using a sample of Ha emitting galaxies. We 
conclude: 

- We confirm only 5 out of 9 emission lines found by 
McCarthy et al. ( 1999). Of the 5 remaining sources, one is a 
broad line AGN, one probably a narrow line AGN. Our conclu- 
sions are therefore based one a very small number of sources. 

- Without reliable inclinations and extinctions, an ensem- 
ble averaged simultaneous inclination and extinction correction 
can be made assuming the relation between extinction and in- 
clination follows that of local galaxies, which is in general a 
small correction because galaxies move approximately along 
the TFR. 

- The most difficult challenge for high redshift TFR studies, 
are the velocity measurements. Besides the observational chal- 
lenges and limitations, there is also the fundamental question if 
what we observe is comparable to what we observe in the local 
galaxies (for example the question if the emission lines extent 
to the flat part of the RC). For the current data set, this leaves 
considerably uncertainties in the measured line widths, which 
cannot be quantified. 

- The high mass end of the local TFR is not sensitive to the 
sample used, therefore the choice of the local reference sample 
is not important for the high mass end of the high redshift TFR. 

- Extending the study of the high redshift TFR to low mass 
galaxies will be very difficult due to the selection effects in ve- 
locity and magnitude: it would require very high spectral reso- 
lution observations of very faint galaxies. 

- Star formation increases the luminosity of galaxies, but 
the effect on log(L) is negligible for high mass galaxies with 
significant older population. Therefore star forming galaxies 
can be TFR-representative for all galaxies, and whether they 
are or not can be checked. 

For the NICMOS galaxies, we measured a ~ 2 mag off- 
set of the z ~ 1 -5 rest frame B starburst TFR with respect to 
the local TFR. This offset is robust against the effects of incli- 
nation and extinction. The linewidths indicate sufficiently mas- 
sive galaxies for TFR studies, the results is therefore also robust 
for the choice of local comparison sample. However, we can- 
not prove or disprove that the all linewidths are due to rotation. 
Moreover, we cannot answer the question if the star burst lumi- 
nosity dominates the total luminosity or not. This is due to the 
extreme nature of the galaxies and the sparseness of data for 
this sample. 

Acknowledgements. We thank the ESO-Paranal staff for executing the 
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